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Abstract 


Initial studies of the sensitivities of estimates of particle fluence, absorbed 
dose j and dose equivalent to fragmentation parameter variations are under- 
taken by using the Langley Research Center galactic cosmic ray transport 
code (HZETRN). The new results, presented as a function of aluminum 
shield thickness , include upper and lower bounds on dose/dose equivalent 
corresponding to the physically realistic extremes of the fragmentation pro- 
cess and the percentage of variation of the dose/dose equivalent as a function 
of fragmentation parameter variation. 


Introduction 

Current methods for estimating interplanetary 
crew radiation exposures resulting from galactic cos- 
mic rays (GCR’s) and their concomitant shielding 
and dosimetry requirements rely on radiation trans- 
port computer codes that describe the interactions 
and propagation of these radiation fields within and 
through bulk matter (refs. 1 through 9). The dose 
and dose equivalent are slowly decreasing functions of 
increasing shield thickness (refs. 3, 4, and 9) because 
of secondary particle production processes whereby 
the heavier GCR nuclei are broken up into nucleons 
and lighter nuclear fragments by nuclear and coulom- 
bic interactions with the shield material. 

As has been discussed elsewhere (ref. 5), the slow 
decrease in dose and dose equivalent with increas- 
ing shield thickness means that relatively small vari- 
ations in predicted doses arising from inaccuracies in 
nuclear fragmentation models may yield large varia- 
tions in estimated shield thickness. Because of the 
paucity of experimental measurements of fragmenta- 
tion cross sections for GCR ions colliding with space- 
craft materials, actual variations in estimates of ab- 
sorbed dose and dose equivalent resulting from the 
use of particular nuclear fragmentation models can- 
not be clearly established at this time. Nevertheless, 
it is possible to use current fragmentation models and 
transport codes to study the sensitivities of dose and 
dose equivalent to variations in these fragmentation 
cross sections. Careful analyses of the extremes of 
the fragmentation event also permit the estimation 
of upper and lower bounds on dose/dose equivalent 
as a function of shield thickness. 

Although currently recommended limits are based 
upon dose equivalent, it is not a good predictor of bi- 
ological risk for heavy ion exposures. The engineer- 
ing design of future, deep-space vehicles will be based 
upon risk models that are more closely related to the 
linear energy transfer (LET) spectra of the radiation 
deposited (ref. 7). Therefore, the effects of variations 


in the fragmentation cross sections on the predicted 
LET spectra will ultimately be of great importance. 

In this work, initial studies of the sensitivity 
of predicted fluences, dose, and dose equivalent to 
variations of the fragmentation cross-section inputs 
into the GCR transport code HZETRN (ref. 8) are 
undertaken. The method used is to simultaneously 
increase or decrease all magnitudes of the fragmenta- 
tion cross-section inputs by the same fixed percent- 
age. At present we have chosen to use 0.50, 0.75, 
1.25, and 1.50 times the nominal values. We have 
also estimated upper and lower bounds on the dose, 
dose equivalent, and LET distributions by consider- 
ing the extremes of the fragmentation event. 

GCR Risk Assessment 

The propagation of these GCR radiation fields 
in shielding is described by a Boltzmann equation 
derived from mass and energy conservation. Its so- 
lutions give particle fluxes and energies everywhere 
within and exiting the boundaries of the target 
medium. For GCR particles, the typically large ion 
kinetic energies allow one to neglect changes in par- 
ticle direction because of collisions. This process 
is called a straight ahead approximation. The one- 
dimensional Boltzmann equation is then written as 
(refs. 2, 8, and 9) 

<t >,( x . E ) = Zj ^j( X ' E ) 

( 1 ) 

In equation (1), 4>i is the flux of type i ions at po- 
sition x with motion along the x-axis and energy E 
(in units of A MeV), <j t is the corresponding macro- 
scopic nuclear absorption cross section (in units of 
cm -1 ), Si is the change in E per unit distance (e.g., 
the stopping power per unit projectile mass), and 
is the cross section (in units of cm -1 ) for producing 
ion i from a collision by ion j. For the GCR trans- 
port problem, the solution is achieved by using the 




first two terms of a perturbation expansion in an it- 
erative procedure which effectively sums subsequent 
generations of reaction products to all orders. The 
result is (refs. 2 and 8) 
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e is a dummy variable, and r'j, the range scale 
parameter relation, is given by 


u i R i {E)=v j Rj{E) (5) 

with Ri{E) denoting the range of ion i with energy E 
and Aj denoting its mass number. The first term on 
the right side of equation (2) represents attenuation 
of the primary ions. The second term on the right 
side represents production and attenuation of secon- 
daries. The distance h is usually chosen to be large 
enough to minimize the number of iteration steps 
while concomitantly minimizing the error resulting 
from the neglect of tertiary production in h. Equa- 
tion (2) can be used to propagate the solution a dis- 
tance h beyond x. Taking the initial point on the 
boundary then allows propagation to any arbitrary 
depth within the interior. Details of the input stop- 
ping powers, nuclear absorption, and fragmentation 
cross sections are given elsewhere (refs. 7 and 8). 

Once the particle fluxes for each GCR ion have 
been determined using equation (2), the absorbed 
dose is computed from 

f oo 

Di(x, >E) = At / Si(E') J>i(x , E?) dE' (6) 

JE 

For risk assessment, the dose equivalent is obtained 
from 


H t (x,> E) = 



S l {E') <t>i{x, E 1 ) dE' 

( 7 ) 


where Qp denotes the quality factor that relates 
absorbed dose to risk. These are LET-dependent 
quantities obtained from ICRP-26 (ref. 10). 


Results 

Sensitivity Studies 

Testing of the sensitivities of the predicted dose/ 
dose equivalent to variations of the fragmentation 
cross sections (<r^ in eqs. (1) and (2)) was accom- 
plished by replacing a l j in equations (1) and (2) by 
pa l j , where the parameter p was varied as desired. 
For this study, we chose the values p = 0.50, 0.75, 
1.00, 1.25, and 1.50. Note that this procedure re- 
sults in a simultaneous, fixed-percentage increase or 
decrease in the magnitudes of all fragmentation cross- 
section inputs. Because of the paucity of experi- 
mental fragmentation cross-section measurements, a 
more realistic approach involving variations of indi- 
vidual cross sections was not considered at this time. 
Such studies should be performed when the experi- 
mental data base for a particular target, such as alu- 
minum, is reasonably well known for many of the 
incident GCR ions. 

The results of the sensitivity studies arc dis- 
played in tables I through III for the solar mini- 
mum GCR environment obtained from the model of 
Adams (ref. 11). The percentage of deviations from 
the nominal values (p = 1.00) is displayed in figures 1 
and 2. For the absorbed dose, which is least sensitive 
to cross-section variations, the maximum deviation 
from the nominal value is approximately 3.5 percent 
for a 50-percent variation in the cross sections. The 
dose equivalent appears to be somewhat more sensi- 
tive in that a 50-perccnt variation in the fragmenta- 
tion cross sections yields a maximum variation in the 
total dose equivalent of approximately 15 percent. 

Table III displays estimated integral fluences by 
particle type for two thicknesses of aluminum shield- 
ing spanning the maximum variations in dose and 
dose equivalent in tables I and II. Note that the 
total flux variations are on the order of 5 per- 
cent for a 50-percent variation in the cross sec- 
tions. The largest variations in the table for a 
given particle type (nearly 70 percent) occur for the 
group of ions with charges between Z — 10 and 
Z = 28. Fluencc variations as large as approximately 
50 percent also occur for the group of ions with 
charges between Z = 3 and Z = 9 and for alpha par- 
ticles (Z — 2). Because of the larger quality factors 
associated with these high-LET particles, the greater 
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variation of the dose equivalent is understandable 
when compared with the variation of absorbed dose, 
even though these particles comprise less than 5 per- 
cent of the total fluence at these shield depths. The 
least variation occurs for protons (4 percent or less) 
and neutrons (8 percent or less) which comprise over 
95 percent of the total transmitted fluence at these 


depths. Variations in LET spectra, which (except 
for bounds on dose equivalent versus LET) are not 
included in this work, may be larger than the dose 
equivalent variations because of dominance by the 
high-LET particles. Such studies, which may be 
relevant to effects on biological systems, should be 
undertaken in the future. 


Table I. Annual Absorbed Dose of Solar Minimum Galactic Cosmic Rays as Function of Fragmentation 
Parameter Variation and Aluminum Shield Thickness 




Annual absorbed dose, cGy, for fragmentation parameters — 


Shield thickness, 
g/cm 2 

0.5 (Tij 

0.75aij 

1 .0 &ij 

1 .25<r t j 

1.50cr ij 

0 

17.06 

17.06 

17.06 

17.06 

17.06 

1 

15.51 

15.47 

15.43 

15.39 

15.35 

2 

15.82 

15.75 

15.68 

15.60 

15.53 

3 

16.01 

15.91 

15.81 

15.71 

15.61 

4 

16.13 

16.00 

15.88 

15.76 

15.64 

5 

16.21 

16.06 

15.91 

15.77 

15.64 

10 

16.26 

16.04 

15.84 

15.64 

15.46 

15 

16.04 

15.78 

15.55 

15.35 

15.16 

20 

15.66 

15.40 

15.16 

14.96 

14.78 

30 

14.68 

14.42 

14.20 

14.02 

13.88 

40 

13.56 

13.31 

13.11 

12.97 

12.85 

50 

12.41 

12.16 

11.99 

11.86 

11.76 


Table IT. Annual Dose Equivalent of Solar Minimum Galactic Cosmic Rays as Function of Fragmentation 
Parameter Variation and Aluminum Shield Thickness 




Annual dose equivalent, cSv, for fragmentation parameters 


Shield thickness, 
g/cm 2 

0 -5(J;j 

0.75crtj 


1.25er,j 

1.50/7^ 

0 

121.7 

121.7 

121.7 

121.7 

121.7 

1 

84.8 

84.3 

83.8 

83.3 

82.9 

2 

82.1 

81.2 

80.3 

79.4 

78.5 

3 

79.3 

78.0 

76.8 

75.6 

74.4 

4 

76.7 

75.1 

73.6 

72.1 

70.7 

5 

74.3 

72.4 

70.7 

68.9 

67.3 

10 

64.8 

62.1 

59.5 

57.2 

55.0 

15 

58.1 

55.0 

52.2 

49.6 

47.4 

20 

53.0 

49.7 

46.9 

44.4 

42.3 

30 

45.4 

42.2 

39.6 

37.5 

35.9 

40 

39.7 

36.7 

34.5 

32.9 

31.7 

50 

34.9 

32.3 

30.6 

29.3 

28.5 
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Table III. Annual Integral Fluences of Solar Minimum Galactic Cosmic Rays as Function of Fragmentation 
Parameter Variation for Representative Aluminum ^Shield Thicknesses 


Particle 

type 

Annual integral fluences, particles/cm 2 , for fragmentation parameters — 

0.5 &ij 

0.75 (Tij 

l.OcJij 

1.25<j t j 

1.5 <Tij 

Aluminum thickness of 20 g/cm 2 

Z — 0 neutron 

1.24 x 10 8 

1.28 x 10 s 

1.32 x 10 8 

1.36 x 10 8 

1.39 x 10 8 

Z — 1 proton 

1.31 x 10 8 

1.33 x 10 8 

1.35 x 10® 

1.36 x 10 8 

1.38 x 10 8 

Z — 2 alphas 

9.26 x 10 6 

8.43 x 10 6 

7.67 x 10 6 

6.99 x 10 6 

6.36 x 10 6 

Z = 3-9 

6.85 x 10 5 

6.24 x 10 5 

5.68 x 10 5 

5.17 x 10 s 

4.70 x 10 5 

Z = 10-28 

1.42 x 10 5 

1.25 x 10 5 

1.09 x 10 5 

9.62 x 10 5 

8.44 x 10 4 

Total 

2.65 x 10 8 

2.70 x 10 8 

2.75 x 10 8 

2.80 x 10 8 

2.84 x 10 8 

Aluminum thickness of 40 g/cm 2 

Z — 0 neutron 

1.57 x 10 s 

1.64 x 10 8 

1.70 x I0 8 

1.74 x 10 8 

1.78 x 10 8 

Z = 1 proton 

1.11 x 10 8 

1.14 x 10 8 

1.16 x 10 8 

1.18 x 10 8 

1.19 x 10 s 

Z = 2 alphas 

7.00 x 10 6 

5.80 x 10 6 

4.80 x 10 6 

3.97 x 10 6 

3.29 x 10 6 

Z = 3 9 

4.69 x 10 5 

3.88 x 10 5 

3.21 x 10 5 

2.64 x 10 5 

2.17 x 10 5 

Z = 10-28 

8.16 x 10 4 

6.30 x 10 4 

4.86 x 10 4 

3.74 x 10 4 

2.88 x 10 4 

Total 

2.76 x 10 8 

2.84 x 10 8 

2.91 x 10 8 

2.96 x 10 8 

3.01 x 10 8 




Figure 1. Percentage of difference from nominal values 
(l.Oaij) of absorbed dose as function of aluminum shield 
thickness. Values are displayed for four different sets of 
input fragmentation parameters (p< 7 tJ ). 


Figure 2. Percentage of difference from nominal values 
(1.0<r t j) of dose equivalent as function of aluminum shield 
thickness. Values are displayed for four different sets of 
input fragmentation parameters ) . 
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Maximum and Minimum Bounds 

Realistic upper and lower limits on the estimated 
dose and dose equivalent resulting from variations in 
the input fragmentation cross sections can be made 
by considering the extremes of a fragmentation event. 
The maximum fragmentation event occurs when the 
projectile nucleus is completely broken up into its 
constituent nucleons. We call this event the low 
LET case because there are no heavy-charged par- 
ticle fragments produced in the reaction. The other 
extreme, the minimum fragmentation event, occurs 
when every collision removes only a single nucleon 
from the fragmenting projectile nucleus. We call this 
event the high-LET case because every reaction pro- 
duces a heavy-charged particle whose mass number 
is one unit less than the mass number of the frag- 
menting nucleus. 

To implement this fragmentation model in the 
Boltzmann equation (eq. (1)), we write the fragmen- 
tation cross section as 

<7ij = TTL^jCj ( 8 ) 

For the low-LET case, where only neutrons and 
protons are produced, 

{ Zj (i = 1 proton) 

Aj - Zj (i — 0 neutron) > (9) 

0 0 , 1 ) J 

For the high-LET case, where only a single neutron 
is removed from the fragmenting nucleus, we use 

mij = (Aj - Zj)/Aj (10) 

for i — j and z* = 0; but for all other events, 
= 0. For the high-LET case where a single 

proton is removed, equation (10) becomes 

m ij = Zj/Aj (11) 

for z = j — l and i — 1; but for all other events, 
mij = 0. Finally, we note that the values of the 
absorption cross sections oj are known to within 
approximately 10 percent (ref. 12). Therefore, we 
also vary <jj by ±10 percent as part of the bounding 
process. 

Figures 3 and 4 display bounds on the absorbed 
dose and dose equivalent obtained using equations (8) 
through (11). For aluminum shield thicknesses be- 
tween 15 and 20 g/cm 2 , differences between the 
upper and lower bounds on the absorbed do_se__ ex- 
ceed 6 percent. For aluminum shield thicknesses be- 
tween 20 and 50 g/cm 2 , differences between the up- 


per and lower bounds on the dose equivalent exceed 
50 percent. 



j 

Aluminum shield thickness, g/cm 


Figure 3. Maximum dose variation in solar minimum galactic 
cosmic rays resulting from consideration of extremes of 
fragmentation event. High LET refers to minimum frag- 
mentation event, and low LET refers to maximum event. 



Figure 4. Maximum dose equivalent variation in solar mini- 
mum galactic cosmic rays resulting from consideration of 
bounds of fragmentation events. High LET refers to min- 
imum fragmentation event, and low LET refers to maxi- 
mum event. 

Bounds on the LET distributions behind 5, 10, 
20, and 30 g/cm 2 of aluminum shielding are shown 
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in figures 5, 6, 7, and 8, respectively. For LET 
values greater than 10 keV/^m, differences exceed 
20 percent behind 5 g/cm 2 of aluminum shielding. 
Increasing the shielding thickness to 10 g/cm 2 of 
aluminum results in LET distributions that differ 
by 50 percent or more. These differences increase 
to greater than a factor of 2 for shield thicknesses 
exceeding 20 g/cm 2 of aluminum. 



Figure 5. Maximum LET distribution variation in solar min- 
imum galactic cosmic rays behind 5 g/cm 2 of aluminum 
shielding resulting from consideration of bounds of frag- 
mentation event. High LET refers to minimum fragmen- 
tation event, and low LET refers to maximum event. 



Figure 6. Maximum LET distribution variation in solar min- 
imum galactic cosmic rays behind 10 g/cm 2 of aluminum 
shielding resulting frorn_ consideration of bounds of frag- 
mentation event. High LET refers to minimum fragmen- 
tation event, and low LET refers to maximum event. 



Figure 7. Maximum LET distribution variation in solar min- 
imum galactic cosmic rays behind 20 g/cm 2 of aluminum 
shielding resulting from consideration of bounds of frag- 
mentation event. High LET refers to minimum fragmen- 
tation event, and low LET refers to maximum event. 



Figure 8. Maximum LET distribution variation in solar min- 
imum galactic cosmic rays behind 30 g/cm 2 of aluminum 
shielding resulting from consideration of bounds of frag- 
mentation event. High LET refers to minimum fragmen- 
tation event, and low LET refers to maximum event. 

Recommended Studies 

The bounds and sensitivity studies reported 
herein comprise an initial investigation of the effects 
of fragmentation parameter uncertainties on calcula- 
tions of galactic cosmic ray exposure. They are by no 
means exhaustive in nature. The quantitative results 
obtained from these studies may be peculiar to the 
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present transport code and its associated inputs. No 
claim is made that alternative transport codes will 
yield these exact results. Further studies that should 
be considered include the following: 

1. Incorporating body self-shielding distribu- 
tions to estimate actual organ exposure 
uncertainties 

2. Estimating exposure uncertainties for other 
shield materials 

3. Expanding the scope of the studies to in- 
clude the newly proposed quality factors from 
ICRP-60 (ref. 13) 

4. Using other GCR environmental models as 
they become available 

5. Carefully studying the effects of parameter 
variations on LET distributions. 

Finally, whenever complete sets of fragmentation 
cross-section data relevant to GCR transport in 
shield materials become available, actual estimates 
of uncertainties in dosimetric quantities should be 
made. 

Concluding Remarks 

Initial studies of the sensitivities of galactic cos- 
mic ray dose and dose equivalent estimates to varia- 
tions in the input fragmentation cross sections were 
carried out for typical thicknesses of spacecraft alu- 
minum shielding. Increases as large as 15 percent 
were noted in dose equivalent when the input frag- 
mentation cross sections were reduced to half of their 
current values. Variations in absorbed dose were 
smaller in magnitude and reflected the predominance 
of high-energy protons in the particle fluence spectra. 
For particles with charge number Z > 2, large vari- 
ations in their fluences were noted. These are im- 
portant because they are the major components of 
particles of high linear energy transfer (LET). 

Upper and lower bounds on dosimetric quan- 
tities were established by considering the physical 
extremes of the fragmentation process. Differences 
between the upper and lower bounds on dose equiv- 
alent exceeded 50 percent for aluminum shield thick- 
nesses greater than 20 g/cm 2 . Above 10 keV/yim, 
differences exceeding a factor of 2 were noted for 
the bounds on the LET distributions. Finally, sug- 
gestions were made in regard to additional studies 
needed on fragmentation parameter sensitivity. 
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